ContourPlot[Cos[2 Pi x] Cos[2 Pi y] == 2 - Exp[-(x - 1)^2/(0.01)^2] Exp[-(y - 1)^2/(0.01)^2], {x, 0, 2}, {y,0, 2}]